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. We study singularity formation in spherically symmetric solitons of 

the charge one sector of the (2+1) dimensional S 2 sigma model, also 
known as CP 1 wave maps, in the adiabatic limit. These equations 
are non-integrable, and so studies are performed numerically on radi- 
ally symmetric solutions using an iterative finite differencing scheme. 
Analytic estimates are made by using an effective Lagrangian cutoff 
outside a ball of fixed radius. We show the geodesic approximation 
^ , is valid when the cutoff is applied, with the cutoff approaching in- 

finity linearly as the reciprocal of the initial velocity. Additionally a 
characterization of the shape of a time slice f(r,T) with T fixed is 
provided. 
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1 Introduction 

In this paper we study a hyperbolic partial differential equation that develops 
a singularity in finite time. 
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The two-dimensional S 2 sigma model has been studied extensively over 
the past few years in ||, 0,0, @, [|J. It is a good toy model for studying 
two-dimensional analogues of elementary particles in the framework of clas- 
sical field theory. Elementary particles are described by classical extended 
solutions of this model, called solitons. This model is extended to (2+1) 
dimensions. The previous solitons are static or time-independent solutions, 
and then the dynamics of these solitons are studied. Since this model is 
not integrable in (2+1) dimensions studies are performed numerically, and 
analytic estimates are made by cutting off the model outside a radius R. 

The model can also be regarded as the continuum limit of an array of 
Heisenberg ferromagnets. 

The S 2 sigma model displays both slow blow up and fast blow up. In 
slow blow up, all relevant speeds go to zero as the singularity is approached. 
In fast blow up the relevant speeds do not go to zero as the singularity is 
approached . The charge 1 sector of the S 2 sigma model exhibits logarithmic 
slow blow up, whereas the charge 2 sector and the similar Yang Mills (4+1) 
dimensional model, both investigated in [|J and [Q, exhibit fast blow up. 

The static Lagrangian density for the S 2 sigma model is given by 

L = J |V0| 2 , 

where is a unit vector field. 

In the dynamic version of this problem, where <fi : IR 2+1 — > S 2 , the La- 
grangian is 

L = I \d t $\ 2 - |V0| 2 . 
Jm. 2 

Identifying S 2 = CP 1 = C U {oo} we can rewrite this in terms of a 
complex scalar field u: 



\d t u\ 2 


|V«| 2 


(1 + \u\ 2 ) 2 


(i + M 2 ) 2 



(1) 



The calculus of variations on this Lagrangian in conjunction with inte- 
gration by parts yields the following equation of motion for the CP 1 model: 

(1 + \u\ 2 )(d 2 u - d 2 x u - d 2 u) = 2u(\d t u\ 2 - \d x u\ 2 - \d y u\ 2 ) (2) 

Here u represents the complex conjugate of u. 



2 



The first thing to identify in this problem are the static solutions deter- 
mined by equation (0). These are outlined in among others. The entire 
space of static solutions can be broken into finite dimensional manifolds M n 
consisting of the harmonic maps of degree n. If n is a positive integer, then 
M n consists of the set of all rational functions of z = x + iy of degree n. For 
this chapter, we restrict our attention to Mi, the charge one sector, on which 
all static solutions have the form 

u = a + (3(z + 7)- 1 . (3) 

In order to simplify, consider only solutions of the form 

pz- 1 

The geodesic approximation says that for slow velocities solutions should 
evolve close to 

z 

Instead of /3(t), we look for a real radially symmetric function f(r,t), and 
find the evolution of: 

f(r,t) 

z 

The differences between the evolution of f(0,t) and that predicted in the 
geodesic approximation, and the deviation of f(r,T) with T fixed has from 
a horizontal line gives us a means to gauge how good the geodesic approxi- 
mation is. 

It is straightforward to calculate the evolution equation for f(r,t). It is: 

The static solutions for f(r,t) are the horizontal lines f(r,t) = c. Here 
c = length scale. In the adiabatic limit motion under small velocities should 
progress from line to line, i.e f(r,t) = c(t). f(r,t) = is a singularity of 
this system, where the instantons are not well defined. We use this to form a 
numerical approximation to the adiabatic limit to observe progression from 
f(r, 0) = Co > towards this singularity. 
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2 Numerics for the CP 1 charge 1 sector model 

A finite difference method is used to compute the evolution of numerically. 
Centered differences are used consistently except for 

dlf + (5) 

r 

In order to avoid serious instabilities in (^) this is modeled in a special way. 
Let 

Lf = r-%r%f = d?f+ 3 -^. 

r 

This operator has negative real spectrum, hence it is stable. The naive central 
differencing scheme on (||) results in unbounded growth at the origin, but the 
natural differencing scheme for this operator does not. It is 
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Questions arise from || about the stability of the solitons for this equa- 
tion. They found that a soliton would shrink without perturbation from 
what should be the resting state. This seems to be a function of the numeri- 
cal scheme used for those experiments. The numerical scheme used here has 
no such stability problems. Experiments show that a stationary solution is 
indeed stationary unless perturbed by the addition of an initial velocity. 

For a full analysis of the stability of this numerical scheme see ||. 

With the differencing explained, we want to derive f(r, t + At). We 
always have an initial guess at f(r, t + At). In the first time step it is 
f(r,t + At) = f(r, t) + t> At with v the initial velocity given in the problem. 
On subsequent time steps f(r, t + At) = 2f(r,t) — f(r,t — At). This can 
be used to compute d t f(r,t) on the right hand side of (|). Then solve for a 
new and improved f(r,t + At) in the differencing for the second derivative 
dff(r, t) and iterate this procedure several times to get increasingly accurate 
values of f(r,t). 

There remains the question of boundary conditions. At the origin f(r,t) 
is presumed to be an even function, and this gives 

/(0,t) = ~/(Ar,t)-^/(2Ar,t). 
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At the r = R boundary we presume that the function is horizontal so 
f(R,t) = f(R-Ar,t). 



3 Predictions of the Geodesic Approximation 

In the time evolution of the shrinking of solitons was studied. They 
arbitrarily cut off the Lagrangian outside of a ball of radius R, to prevent 
logarithmic divergence of the integral for the kinetic energy, and then analyze 
what happens in the R —>■ oo limit. In || the problem of the logarithmic 
divergence in the kinetic energy integral is solved by investigating the model 
on the sphere S 2 . The radius of the sphere determines a parameter for the 
size analogous to the parameter R for the size of the ball the Lagrangian is 
evaluated on in ||. 

Our calculations here follow those done in ||. 

Equation (HD gives us the Lagrangian for the general version of this prob- 
lem. In the geodesic approximation or adiabatic limit we have 

(3 

u = — 

z 

for our evolution. If we restrict the Lagrangian to this space we get an 
effective Lagrangian. The integral of the spatial derivatives of u gives a 
constant, the Bogomol'nyi bound, and hence can be ignored. If one integrates 
the kinetic term over the entire plane, one sees it diverges logarithmically, so 
if (3 is a function of time, the soliton has infinite energy. 

Nonetheless, this is what we wish to investigate. We cannot address the 
entire plane in our numerical procedure either, hence we presume that the 
evolution takes place in a ball around the origin of size R. If (5 — f(r,t) 
shrinks to in time T, we need R > T. Under these assumptions, up to a 
multiplicative constant, the effective or cutoff Lagrangian becomes 

f R r%p 

Jo FT7¥ 

which integrates to 



Since the potential energy is constant, so is the purely kinetic Lagrangian, 
and 



r»2 

In I 1 + 4 



R 1 



c 

2"' 



P) P + R\ 
with c (and hence c 2 /2) a constant. Solving for <9 t / we obtain 

dtf = 



d2 

In I 1 + 5r 



i? 2 



(6) 



Since we are starting at some value fo and evolving toward the singularity 
at / = this gives: 



f(o,t) 




(7) 



The integral on the right gives ct. The integral on the left can be evaluated 
numerically for given values of R, f and /(0, t). A plot can then be generated 
for ct vs. /(0, t). What we really are concerned with is /(0, t) vs. t, but once 
the value of c is determined this can be easily obtained. One such plot with 
fo = 1.0, R = 100 of /(0, t) vs ct is given in Figure [L| This curve is not quite 
linear, as seen by comparison with the best fit line to this data which is also 
plotted in Figure |I| The best fit line is obtained by a least squares method. 
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4 Evolution of the Origin Results 



The computer model was run under the condition that f(r, 0) = fo with 
various small velocities. The initial velocity is d t f(r, 0) = v o, other input 
parameters are R = r max , Ar and At. 
Evolution of f(0,t) 

The primary concern with the evolution of the horizontal line is the way 
in which the singularity at f(0,t) = is approached, because once again as 
r — > oo equation (||) reduces to the linear wave equation 



and so we expect the interesting behavior to occur near r = 0. The model 
is run with initial conditions that f(r, 0) = fo, d t f(r,0) = vq. Other input 
parameters are R = r max , Ar and At. 

The evolution of the initial horizontal line seems to remain largely flat 
and horizontal, although there is some slope downward as time increases. 
This is shown in figure [|. 

We track f(0,t) as it heads toward this singularity, and find that its 
trajectory is not quite linear, as seen in figure |3|. This is suggestive of the 
result obtained in the predictions for this model. 

We want to check the legitimacy of the result from equation (|7]) in chapter 
3.2. This requires a determination of the parameters R, the cutoff, and c. 
We already have fo and f(0,t). To determine R and c, observe from equation 



d 2 j = a?/, 



d) that: 








dtP 



8 



The plot of ln(/) = ln(/(0,t)) vs l/d t f 2 = l/d t f(0,t) should be linear 
with the slope m = 2/c 2 and the intercept b = (21n(i?) — l)/c 2 . Such a 
plot is easily obtained from the model, and given slope and intercept, the 
parameters c and R are easily obtained. 

Figure § is a plot of ln(/(0, i)) vs. l/<9 t / 2 (0, t), with initial conditions 
Ar = 0.01, At = 0.001, f = 1.0 and v = -0.01. It is easily seen that 
although the plot of ln(/(0, t)) vs. l/d t f 2 (0,t) is nearly straight, it is not 
quite a straight line. This may indicate that the values of R and c are 
changing with time. 

The best fit line y = mx + b has slope m = —2810 and b = 10200. We 
have 

V m 

and 

R = exp ( + 1) . 

This gives values c = 0.0267 and R = 62.1. Using these values of c and R 
in the calculation of equation (0), we obtain the plot of f(0,t) vs t given 
in Figure [|. This is overlayed with the model data for f(0,t) vs. t for 
comparison. These two are virtually identical. 

This shows that the phenomenon of cutting off the Lagrangian outside 
of a ball of radius R is not just an artifact of necessity because the full 
Lagrangian is divergent, but an inherent feature of this system. 

Table [I] contains the data for c and R vs. change in the initial velocity 
Vq, under the initial conditions /o = 1.0, Ar = 0.01 and At = 0.001. The 
data for R varying with 1/vq fits well to the line 

y = 0.5407a; + 6.032. 

This fit is shown in figure [| A linear fit makes sense, since as the velocity 
tends toward zero, we expect the cutoff to head toward infinity. 

Table H containing the data for c and R vs. change in f , the initial 
height. The parameter R varies close to linearly with / , while the parameter 
c remains nearly constant, changing by less than 7% over the course of the 
runs. 
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Time Profiles of Data 
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Figure 2: Time slices f(r,T) 



Table 1: c and R vs. v . f = 1.0. 



^0 


c 


R 


-0.005 


0.0145 


115 


-0.00667 


0.0187 


89 


-0.01 


0.0263 


53 


-0.0133 


0.0342 


49 


-0.02 


0.0485 


34 


-0.03 


0.0683 


25 


-0.05 


0.104 


17 


-0.06 


0.121 


15 
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Evolution of Origin 




50 60 70 80 90 100 110 120 

Time in seconds 



Figure 3: Evolution of f(0,i) is not quite linear. 



Table 2: c and R vs. /q. v = —0.01. 



fo 


c 


R 


1.0 


0.0267 


62 


2.0 


0.0263 


108 


3.0 


0.0260 


150 


4.0 


0.0259 


190 
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Figure 5: Predicted course of /(0, t) from equation (0) and actual course 
/(0,f)vs. t. 
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R vs 1/vO and best linear fit to the data 
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Figure 6: R vs 1/vO and the best fit line. 
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Characterization of time slices f(r,T) 

Making a closer inspection of the time profiles f(r,T) with T fixed as 
in Figure 0, one may observe that the initial part of the data is close to a 
hyperbola as seen in Figure [7]. The best hyperbolic fit is determined by a 
least squares method. 

The equation for the hyperbola is 



b 2 a 2 

One would naturally ask about the evolution of the hyperbolic parameters 
a and b with time, however, neither of these is particularly edifying. A simple 
calculation shows that k should follow /(0, t) closely if b is small, as it is. 

The evolution of — b/a gives the slope of the asymptotic line to the hy- 
perbola, and this evolution is close to linear as seen in Figure |[ 
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Hyperbolic parameter -b/a as a function of time with best fit line 




gure 8: Plot of hyperbolic parameter — b/a vs time, f = 1.0, v = — 0. 
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5 Conclusions 

In ||, solitons are found to be numerically unstable. A solution of the form 

t 

z 

shrinks spontaneously under their numerical procedure. This does not occur 
in our numerical implementation of the S 2 sigma model. The static solutions 
do not evolve in time unless given an initial rate of shrinking. Further, stabil- 
ity and convergence analysis of two of the numerical procedures is provided 
in§- 

We use a method analogous to that in Q of cutting off the Lagrangian 
outside of a ball of radius R, and we find an explicit integral for the shrinking 
of the soliton, dependent on two parameters: c which is a function of the 
kinetic energy, and R which is the size of the ball on which we evaluate 
the Lagrangian. Once these are specified, this integral gives the theoretical 
trajectory of the soliton. We find an explicit integral for the shrinking of a 
soliton where the Lagrangian is cut off outside of the ball of radius R, and 
when R is calculated (along with the parameter c), the shrinking seen in the 
numerical model matches that predicted. The dependence of R on the initial 
conditions appears to be linear with the initial velocity. The validity of the 
geodesic approximation with a cutoff R is shown by this, and the cutoff is 
an essential part of this system. 

In addition to this, the shape of a time slice f(r,T), with T fixed, is 
characterized by hyperbolic bumps at the origin. 
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